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Abstract 

A new spectrum generating algebra for a unified description of rotations and 
vibrations in polyatomic molecules is introduced. An application to nonlinear 
X 3 molecules shows that this model (i) incorporates exactly the relevant point 
group, (ii) provides a complete classification of oblate top states, and (iii) treats 
properly both degenerate and nondegenerate vibrations. 
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1 Introduction 


The progress currently witnessed in experimental techniques for the spectroscopy of 
large molecules motivates the development of new theoretical tools to interpret and 
guide such measurements. In view of the proliferation of parameters needed in Dunham 
expansions and the difficulties encountered in solving a Schrodinger equation with 
interatomic potentials for polyatomic molecules, it has been suggested to use algebraic 
methods in which the hamiltonian is expressed in terms of elements of a Lie algebra. 
The key ingredient is the choice of a suitable spectrum generating algebra. 

At present there exist two versions of this method. In the first version ill , rota¬ 
tions and vibrations are treated simultaneously, and hence rotation-vibration couplings 
are built in from the outset 0 . In this approach, a molecule with n atoms is described 
in terms of n — 1 coupled U (4) groups, one for each independent relative coordinate. 
Each Ui( 4) algebra is realized in terms of a set of four vibron operators (< 7 *, 7 q): the 
three components of a dipole (or 7T-) boson with L p = 1“ and a scalar (or a-) bo¬ 
son with L p = 0 + , i.e. one cr-boson for every 7 r-boson. The scalar boson does not 
represent an independent degree of freedom, but is merely introduced to conveniently 
handle anharmonicities by compactifying the model space, i. e. the number of bosons 
N t = n ai + n ni , is conserved for each i — 1,..., n — 1 separately. This version has been 
applied successfully to polyatomic molecules with n — 2,3,4 atoms, but it encounters 
difficulties in describing bent molecules with degenerate vibrations, such as for example 
in nonlinear X 3 molecules p. ff|. 

In the second version m 0 rotations are ignored, while vibrations are treated in 
terms of coupled one-dimensional anharmonic oscillators. A separate set of coupled 
17(2) groups is introduced for each type of vibration (stretching and bending). The 
use of symmetry adapted operators ensures the correct transformation properties of 
the eigenstates under the relevant point group 0. Since there is no explicit coupling 
between rotations and vibrations, this scheme becomes particularly useful for molecules, 
in which the rot at ion-vibration coupling is negligible. I 11 this approach each type 
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of vibration is treated separately with different interactions, although for degenerate 
vibrations, such as for example in symmetric top X 3 molecules, stretching and bending 
vibrations can belong to the same irreducible representation (irrep) of the relevant 
point group. 

In this note, we propose an alternative scheme for the description of n-atomic 
molecules, in which the relevant point group symmetry is taken into account exactly 
and which contains all vibrations (stretching and bending) and rotations in a single 
algebraic framework. We introduce a dipole boson for each independent relative co¬ 
ordinate, and a single scalar boson. This leads to a spectrum generating algebra of 
U(k + 1), where k = 3 (n — 1) is the total number of rotational and vibrational degrees 
of freedom. By construction, only the total number of bosons, N = n a + n n with 
n n = n iu is conserved. For diatomic molecules (n = 2) we recover the U (4) vibron 
model (![]. For triatomic molecules (n = 3) we obtain a U( 7) model whose building 
blocks are a scalar boson and two dipole bosons (a, 7 Ti, 7 t 2 ). We present this scheme 
by studying symmetric (oblate) top X 3 molecules which form the simplest nontrivial 
example of a bent polyatomic molecule with a degenerate vibration. We compare our 
results with those obtained in a 17(4) x f/(4) model [2j]. 

2 Point group symmetry 

The embedding of discrete point group symmetries in an algebraic model relies on a 
geometric interpretation of the vibrons. For nonlinear rigid molecules it is convenient 
to make use of the established isomorphism between the molecular symmetry group 
and the relevant point group (?]]. The elements of the former consist of permutations 
of identical nuclei with or without inversion (parity). The parity of the bosons is 
well defined and the transformation properties under permutations can be realized in 
terms of finite rotations among the bosons. For bent X 3 molecules with point group 
symmetry D 3 h, we use the isomorphism between the S 3 permutation group and the 
point group D 3 (c D 3h ). The parity and S 3 labels are equivalent to the classification 
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under D 3 h- We associate the two types of dipole bosons, Ti\ and 7t 2 , with the relative 
Jacobi coordinates 


P 


1 

71 


(xi - X 2 ) , 
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(fi + x 2 


2 ^s) , 


( 1 ) 


and their conjugate momenta, which have well defined transformation properties under 
permutations (x, denotes the coordinate of the i-th atom). The elements of the group 
S 3 can be expressed in terms of the transposition P(12) and the cyclic permutation 
P(123) whose matrix representation in the (cd, 7,m; 71 2 ,m) basis is || 
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The cr-boson is a scalar under the permutation group, whereas the two 7r-bosons trans¬ 
form (for each projection m — 0, ±1) as the two components ( M p , M\) of the two- 
dimensional irrep (M) of S 3 . In general, the symmetric (S), antisymmetric (A), and 
mixed symmetry (M) classes of S 3 , can equivalently be labeled by the irreps of the 
isomorphic point group D 3 , as A\, A 2 and E, respectively. 

The transformation properties under S 3 of all operators of interest follow from 
those of the building blocks. In particular, the most general one- and two-body U(7) 
hamiltonian which is a scalar under S 3 as well as rotationally and parity invariant is 
found to be 


H = e a cdcr — (7 • 7Ti + 7 ' ^ 2 ) + Uq cdcdcnr — u\ cd(7 • ffi + 7 ' ^ 2 ) o' 

+ v 0 [(7 ■ 7 + 7 • 7) crcr + cr t cr t ( 7 fi • 7 Ti + 7 r 2 • 7 t 2 ) 

+ J2 C A (77 - 77) (A) ■ (TTl^l - ^27T2) (A) + 4 (77) (A) ■ (7T27Ti) (A) ] 

A=0,2 

+Ci (77) (1) • (7T27Tl) (1) + J2 W * (77 + 77) (A) • (^1 + n 2 7t 2 ) {X) , (3) 

A=0,2 


with 7 iy m = (—and i = 1,2. The corresponding eigenstates are labeled by 
the total number of bosons N and, by construction, have good angular momentum, 
parity and permutation (point group) symmetry. 
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3 Geometry 


The hamiltonian of Eq. (||) is expressed in terms of abstract algebraic interactions. A 
more intuitive geometric visualization can be obtained by using mean-field techniques 
to study the geometric properties of the {7(7) model. For a system of bosons the 
variational wave function takes the form of a coherent state || which is a condensate 
of N bosons, 


with 



(4) 


bl = (1 + R 2 ) cr t + ri 7r| , + r 2 (cos dirlz + sin d7r 23 ,) 


(5) 


where R 2 = rf + r^. The condensate is parametrized in terms of two (real) coordinates, 
r\ and r 2 , and an angle 6, (ri, r 2 > 0 and 0 < 9 < tt). The two vectors ri and r 2 span the 
xz-plane. We have chosen the £-axis along the direction of ri, and r 2 is rotated by an 
angle 6 about the out-of-plane y- axis, r 1 -r 2 = r x r 2 cos 6 . The expectation value of the 
hamiltonian of Eq. ([lj) in the condensate defines a classical energy surface -E(ri,r 2 ,0). 
The equilibrium shape is determined by minimizing the energy surface with respect 
to ri, r 2 and 0. The nonlinear rigid and stable equilibrium shape is characterized by 
r\ = r 2 and 6 = tt/2. These two conditions are precisely those satisfied by the Jacobi 
coordinates of Eq. ([!]) for an equilateral triangular shape, and suggests to associate the 
algebraic shape parameters in Eq. (Q) with these coordinates, i.e. r\ p and r 2 A . 


4 Excitations 


We consider separately the vibrational and rotational excitations described by the 
hamiltonian in Eq. (||). With the techniques introduced in [Jl(Jl for the nuclear inter¬ 
acting boson model, an arbitrary U(7) hamiltonian can be decomposed uniquely into 
intrinsic (vibrational) and collective (rotational and rotation-vibration coupling) parts, 


77 = 77 ; 


int 


77 . 


coll 


( 6 ) 
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4.1 Vibrations 


The intrinsic part by definition annihilates the equilibrium condensate and has the 
same shape for the energy surface as the original hamiltonian. For the rigid triangular 
equilibrium shape, characterized by rq = r 2 and 6 = tt/2 , we find 


^int - 


6 (R 
+6 ( 


2 <7 V — 7r| • 7r| — 4 ' 4 


)(« 2 


<7(7 — 7iq ' 7iq — 7T 2 ■ 7T 2 


4 • 7rj 


4 • 4 


7Ti • 7Ti - 7T 2 • 7T 2 ) +4(4'4) (^2 ' ^l) 


(7) 


To get more insight into the vibrational structure, we perform a normal-mode analysis 
to the above hamiltonian. This is done by introducing a set of orthonormal deformed 
vibrons, consisting of the condensate boson of Eq. (j5|) with rq = r 2 and 6 = n/2, and 
six additional bosons which represent excitations of the condensate. The normal modes 
can be found by rewriting the intrinsic hamiltonian in terms of the deformed vibrons 
and replacing the condensate bosons by their classical mean-field value \/N ||TT| , [T^ . 
As a result we find to leading order in N 


~ e i blpu + e 2 {b\b v + bl,b w ) + 0(1/W) , (8) 

with eigenfrequencies e\ = 4£i-R 2 and e 2 = 4^ 2 A 2 (1 + A 2 ) -1 . This identifies the 
deformed bosons that correspond to the three fundamental vibrations: a symmetric 
stretching (w), an antisymmetric stretching ( v ) and a bending vibration ( w ), 


4 

4 

bl 


(1 + R 2 ) 1/2 [-Ra^ + (4, 2 + 4.x 


(-4, 2 + 4J/4 2 , 
(4,s + 4, 2 )/4 2 . 


(9) 


The first two are radial excitations, whereas the third is an angular mode which cor¬ 
responds to oscillations in the angle 9 between the two Jacobi coordinates (iq, 
z/ 2a and U 2 b, respectively in the usual spectroscopic notation). The angular mode is 
degenerate with the antisymmetric radial mode. This is in agreement with the point 
group classification of the fundamental vibrations for a symmetric X 3 configuration 
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and shows that H j n ^ describes the vibrational excitations of an oblate symmetric 
top. Intrinsic states representing excited vibrations (zq, u l 2 ), are obtained, for large N, 
by replacing the condensate bosons in Eq. (|||) by the appropriate number of deformed 
bosons of Eq. (S|). Accordingly, the intrinsic state for such vibration (having +1 pro¬ 
jection of the vibrational angular momentum on the symmetry y-axis) takes the form 
(^,i) ( " 2 + 0 / 2 (^,-i ) ( " 2_0/2 (^)" 1 (^) 7V_ " 1_ " 2 I 0 ), where b{ ±1 = (b\±ibl)/V2. 

By construction, has an exactly degenerate ground band whose rotational 

members are obtained by projection from the equilibrium condensate. Its excited states 
also tend to cluster into bands. Eq. @ shows that in the large N limit the vibrational 
spectrum is harmonic. For finite values of N the correction terms of order 0(1/\/N) 
give rise to anharmonicities. These can be studied numerically by diagonalizing 
in a convenient basis. 

A similar analysis [|T 2 | of the ^-invariant intrinsic (one- and two-body) hamiltonian 
of the 1/(4) ® 17(4) model yields that, in this case, the two radial modes associated with 
the Jacobi coordinates are uncoupled and have the same frequency, in disagreement 
with the point group classification of the normal vibrations of a symmetric X 3 shape. 
The difference with the U( 7) model can be traced back to the occurrence of terms in 
the hamiltonian of Eq. (0) in which only the total number of dipole bosons is conserved. 
These terms cannot appear in a U(4)®U(4) description which requires that each type of 
bosons is conserved separately. A U(2)®U(2) description will yield only the symmetric 
Vi(Ai) and the antisymmetric v 2a (E) stretching vibrations. The second member of the 
degenerate vibration (the v 2 b(E) bending) requires a separate treatment. 


4.2 Rotations 

On top of each vibrational excitation there is a whole series of rotational states. In a 
geometric description the rotational excitations are labeled by the angular momentum 
L and its projection K (= 0,1,...) on the threefold symmetry axis, parity P = (—) A 
and the transformation property t under the point group. In the present algebraic 
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model the rotations (and rotation-vibration couplings) are described by the collective 
part of the hamiltonian. By construction, H QQ jj — H — H- U] ^ consists of interaction 
terms which do not affect the shape of the energy surface ITU . Discarding TV- 
dependent terms that do not contribute to the excitation spectrum, we find 


with 


tfcoll _ Kl ' A\ + A 2 ■ A 2 ) + K'2 (B 1 • B\ + P 2 ■ P 2 ) 

T /?3 L • L + K 4 By ■ Ky , 

M = i ( 7 r\a + cr t 7fi)( 1) , A 2 — i (7r|cr + cr t 7f 2 ) (1) , 

B\ = (4^2 + 7r|7Ti) (1 ) , B 2 = (Trjffi - 7r 2 7r 2 ) (1) , 

L = \/2 (7r|#i + 7T^7f 2 ) (1) , ky = -ik 3 (7r|7f 2 - 7T^7fi) (0) . 


( 10 ) 


( 11 ) 


The angular momentum L commutes with any rotational-invariant hamiltonian. The 
hamiltonian of Eq. (||) commutes also with the operator K y . Consequently the re¬ 
sulting eigenstates \L p ,K y ) have good angular momentum L, parity P and K y . The 
states | L p , ±K y ) are degenerate since P( 12) \L P , K y ) = (—) Ky I L p , —K y ) and P( 12) 
commutes with any S 3 invariant hamiltonian. The operator K y is related to the cyclic 
permutation P(123) = exp (—i2irK y /3) and the absolute value K y determines the 
S 3 ~ D 3 symmetry [Q, 

0 (mod 3) ior Ai,A 2 

1, 2 (mod 3) for E 


\Ky\ = 


( 12 ) 


States with well defined permutation symmetry ( t ) can be formed by taking the two 
linear combinations | L p , \K y \) oc 1±P(12) | L p , K y ), which can be used to distinguish 
between t — A\ and A 2 for \K y \ = 0 (mod 3), and between the components E p and E\ 
for \K y \ = 1,2 (mod 3). 

By examining the projection of intrinsic states on the symmetry axis, it is possible 
to show that for a (zq,z4) vibration, the algebraic quantum number \K y \ and the 
geometric K are related by 


\K y \ = \K±2l\ . 


( 13 ) 







The +(—) sign corresponds to the +/(—/) projections of the vibrational angular mo¬ 
mentum along the symmetry axis. The \K y \ quantum number is analogous to the G 
quantum number 


G = \K±l\ , 


(14) 


considered by Watson Jl5|. Since both the G and \K y \ labels are defined (mod 3), they 
provide an equivalent classification scheme. It is important to note that, unlike for 
l = 0 where K y \ — G — \K\, for l > 0 there are two possible \K y | or G values for each 
K. Thus K y (or G) is an additional quantum number needed to supplement the D% 
and parity labels for a complete classification of the oblate top states. As an example, 
in if-vibrations of the type (u 1 ,u l 2 =1 ) the two L^ levels with K = 3 are distinguished 
by | K y \ = 5 (G = 2) and \K y \ = 1 (G = 4). 


From the above discussion it is evident that the last two terms in Eq. (|T0]) com¬ 
mute with the hamiltonian of Eq. ([3|) and thus correspond to exact symmetries. Their 
eigenvalues K 3 L(L + 1) + k 4 K y , are similar in form to those of a symmetric top. All 
rotational states with K 7 ^ 0 are split in two, with different \K y | assignments in accord 
with Eqs. (12) and (13). The splitting between +1 and —l levels is 8 k 4 /A', i.e. increases 
linearly with K. The K\ and k 2 terms in Eq. (|lOj) do not commute with the intrin¬ 
sic hamiltonian of Eq. (^). Therefore, in addition to shifting and splitting the bands 
generated by H^, they can also mix them and hence contain the rotation-vibration 
couplings. Their effect on the spectrum can be studied numerically. 

The one- and two-body U( 7) hamiltonian presented so far does not exhibit /-type 
doubling and consequently its calculated Ai and A 2 levels occur in degenerate doublets. 
This degeneracy can be lifted with the inclusion of higher-order terms which break the 
K y symmetry. An example of such a three-body term is considered below. A similar 


situation is encountered in Watson’s effective hamiltonian [|T5|, whose major terms are 
diagonal in the quantum number G, and small higher-order correction terms (with 
AG = 6 ) are important for the splitting of the (Ai, A 2 ) doublets. 
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4.3 Applications 


Among the new features in the proposed U( 7) spectrum generating algebra, as com¬ 
pared to the U\ (4) ® 172(4) scheme, is the exact account of the point symmetry. This is 
particularly important for degenerate vibrations present in symmetric top X 3 molecules. 
With that in mind, we apply the present formalism to the vibrational spectrum of H-/. 
This molecule, with a triangular equilibrium shape, has been studied extensively ]16 


and can serve as a test ground for different spectrum generating algebras. It was shown 


l2| that the 1/(4) <8)17(4) model encounters some difficulties in the description of the 
vibrational spectrum of H^" 0 . Our goal here is to examine to what extent these diffi¬ 
culties arise from the lack of an exact treatment of the molecular point group symmetry. 
For that purpose we use a simple 1/(7) hamiltonian 


H = if int (fi, 6, R 2 ) + K A k y -ky + CT^ -T , 


(15) 


where is given in Eq. (7) and = (<7rj -7r| — ^2'^2) 7r 2+2(7r|-7r2)7r|. The parameters 
are determined from a fit to the estimated band origins in hf/ |L7|. In the fit we 
minimize J2i w i [E% (calc)—E)] 2 with weights tty = 1.0 for the two fundamentals (the only 
band origins known from experiment) and wy = 0.1 for the other levels. As shown in 
Table 1, the 1/(7) calculation shows a significant improvement over the U± (4) <8>I/ 2 (4) fit, 
thus highlighting the importance of incorporating the discrete point group symmetry 
in the algebraic description. It is clear, however, that the simple 1/(7) hamiltonian (15) 
used here is not sufficient to reach the accuracy of ab initio calculations. Adding more 
interactions ( e.g . rot at ion-vibration terms and/or additional higher-order terms in the 
algebraic hamiltonian) can improve the spectroscopic accuracy of the fit at the expense 
of more parameters. The need for such terms is expected in view of the known large 
rotation-vibration couplings in Hj/, the non-Born Oppenhcimer corrections in the H^" 
potential [H| and the presence of higher-order terms in ab initio potential surfaces (31 
terms fl7H). 
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5 Summary and conclusions 


In this work we have introduced a U( 7) spectrum generating algebra for the descrip¬ 
tion of triatomic molecules. Particular emphasis was put on the ability of the model to 
adequately treat degenerate vibrations as they occur in oblate top X 3 molecules, by im¬ 
posing the discrete point group symmetry on the hamiltonian. A normal mode analysis 
revealed that the characteristic pattern of fundamental vibrations of a symmetric bent 
X 3 molecule is recovered in this model. Algebraic terms affecting the vibrations, rota¬ 
tions and their coupling have been identified. The U( 7) algebra provides a quantum 
number \K y \ (the analog of the quantum number G ), needed for a complete classifi¬ 
cation of oblate top states. Although we have focussed the discussion to oblate top 
X 3 molecules, the U( 7) model can be applied to any triatomic molecule, both linear 
and bent. In each case the results of a normal mode analysis are in agreement with 
those dictated by the relevant point group All these features of the U( 7) model are 
necessary ingredients for a proper description of rotations and vibrations in triatomic 
molecules. 

We have compared the present model with two other algebraic approaches. The 
main difference between the models lies in the choice of the model space. In the U( 7) 
model the distribution of quanta among the two dipole degrees of freedom is determined 
dynamically by the hamiltonian, whereas in the 17(4) ® 17(4) model (or its simplified 
17(2) <8) 17(2) version) the number of bosons in each mode is restricted separately. The 
existence of generators in the 17(7) algebra that mix the two types of dipole bosons, 
is the new ingredient that enables an exact treatment of the point group symmetry 
and a proper description of degenerate vibrations. Unlike in the 17i(4) ® 17 2 (4) case, 
in the 17(7) scheme there is no need to use non-analytic (absolute value) interaction 
terms to describe bent molecules. On the other hand the 17(7) algebra does not have 
a .SO (4) ® 5*0(4) dynamical symmetry, which led to considerable simplification in 
the 17(4) ® 17(4) description of linear molecules. In a fit to the band origins of we 
showed that a simple U (7) hamiltonian with the proper point group symmetry improves 
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considerably the previous U\(4) g) 172(4) description of triangular X 3 molecules. 

Refinements of the simple [/(7) hamiltonian are needed to achieve spectroscopic 
accuracy. The choice of rotation-vibration terms and of additional higher-order alge¬ 
braic terms may be guided by examining corresponding terms in Watson’s effective 
hamiltonian. The inherent simplicity in algebraic methods is potentially important for 
polyatomic molecules, where are ab initio calculations are difficult to perform. A first 
step in this direction would be to extend the current procedure to larger molecules, 
resulting in a spectrum generating algebra of U (3 n — 2) as a candidate for a spectrum 
generating algebra for vibrations and rotations in n-atomic molecules with prescribed 
point group symmetries. Work along these lines, as well as on the rotational aspects 
of these models, is in progress and will be reported separately. 
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Table 1: band origins in cm 


Oi,r4) 

t 

DMT“> 

U(7) b) 

*7(4) 8) U( 4) c) 

I II 

(0,0°) 

A\ 

0 

0 

0 

0 

(0,! 1 ) 

E 

2521 

2522 

2506 

2506 

(1,0°) 

A\ 

3178 

3184 

3181 

3174 

(0,2°) 

A\ 

4778 

4770 

4928 

4929 

(0, 2 2 ) 

E 

4998 

5007 

4928 

4929 

(l,! 1 ) 

E 

5554 

5543 

5603 

5597 

(2,0°) 

A\ 

6262 

6233 

6233 

6243 

(0.3 1 ) 

E 

7006 

7019 

7266 

7269 

(0,3 3 ) 

A\ 

7283 

7276 

7266 

7269 

(0,3 3 ) 

A 2 

7493 

7492 

7266 

7269 


“•* From 


171 


^ Calculated with N = 40, R 2 = 1.65, £i = 


12.3 cnr 1 , £ 2 = 21.3 cm -1 , k 4 = —14.8 cm -1 
and £ = —0.38 cnr” 1 
c ) From [Q] 
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